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A software package was developed to provide the inputs fora 
passive single-sensor acoustic tracker. This required a computation 
of sound propagation velocity and resolving high-resolution frequency 
information from raw acoustic data. A selected frequency was then 
automatically tracked from an up-doppler to a down-doppler condition. 
This frequency information coupled with associated bearing information 


was then used to provide continuous target fixes. 
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I. INTRODUCTION 


Continuous fixing of an underwater sound source by passive means 
has traditionally been carried out by using multiple sensors. Single 
sensor fixing requires varying power coefficients, frequencies, and/or 
.bearing.changes and usually produces an estimated position with an 
accuracy that varies greatly with the type of sensor and processing 
utilized and the skill of the evaluator. An accurate single sensor track 
would be advantageous in that positional errors could be minimized with 
only one buoy or array required, providing greater flexibility in that a 
series of sensors are not required to be monitored, and a decrease in 
buoy expenditures could be expected. To provide the necessary frequency 
data, a software package was developed that computed the sound propaga- 
tion velocity, detected a discrete frequency signal and automatically 
tracked that signal from a full up-doppler through a full down-doppler 
situation. This frequency output combined with associated bearing in- 
formation, and a value for the sound propagation speed with an estimate 
of the speed of the target provides sufficient input data for single sensor 


fixing. 
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Ir THE TERRERING PROBLEM 


A. STATEMENT OF THE PROBLEM 

Chapter 3 of [1] developes the required input quantities for passive, 
single-sensor fixing. These values are an estimate of target speed, the 
velocity of sound propagation in the medium, the doppler shifted fre- 
quency and a bearing from the sensor to the target. The cstimate of 
target speed would have to be determined by tactical considerations at 
the time and location of signal acquisition, i.c., target transiting or on- 
station, signal strength, detection ranges, etc. The bearing would be 
provided by a standard DIFAR buoy or an array. In this problem the 
standard deviation of bearing measurement noise was varied from 4 Lo 
30%, Methods of providing the sound propagation velocity and the doppler 


shifted frequency were then devised. 


В. DEVELOPING THE INPUTS 

Equation 1.1. of [2] provides a good approximation of sound speed 
In scaswaler as a function of temperature, salinity and pressure. It is 
repeated here as equation 1. 

c(T,S,z)-144944. 62 T-. 054T^ 41. SSE OI (1) 

where c is in meters per second; T is in degrees centigrade, S is in 
parts per thousand, and z is in meters. Converting this equation so 
that c is in yards per second, T is in degrecs Farenheit and depth is in 


fect, yields cquation 2. 


11 





c(T, S, z)=1476. 143. 97T-. O18T“ 41. 42(S.35)+. 0057z (2) 
If the target depth is known, the median depth of the target and hydro- 
phone should be used. Otherwise, use the hydrophone depth. Salinity 
would normally be set to 35ppt unless operating in water of a known 
different salinity. The resultant output from this equation would be set 
equal to the velocity of sound propagation, VP, and used as an input 
value to the target tracker program. 

The frequency requirements called for a high-resolution value 
which dictated a relatively long time window. <A twenty second window 
was used which resulted ina .05Hz FFT resolution. This resolution 
was further enhanced by interpolation using the relative values of the 
power coefficients associated with the discrete frequency being tracked. 
In order to maintain contact on the same discrete frequency line, a 
frequency tracker algorithm was devised. This provided the necessary 


frequency data for the target tracker program. 
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II. SOFTWARE IMPLEMENTATION 


A. SOUND SPEED 

This process was a simple matter of implernenting the equation for. 
sound speed in seawatcr, cquation 2. Providing the chosen depth, ob- 
served temperature, and salinity results in a best estimate of the actual 
sound speed in the operating area. Updating this valuc would be necessary 


as the opcrating area or cnvironmental conditions changed. 


Б. “ЕБКЕОТОЕМСҮ 
1. Parameters for Frequency Computations 

In order to provide realistic frequency inputs, a magnctic tape 
was obtained that contained analog, raw acoustic signals representing 
arabient noise plus an underwater sound cmitter with various discrete 
frequencies. The doppler shift ofa discrete frequency signal is directly 
proportional to the rest frequency of the signal. Thus doppler shift 
measurements are more casily obtained from higher frequencies. 
However, signal attenuation in the medium is much higher at highcr 
frequencies. At present, frequencies in the 150Hz to 300Hz range will 
yicld usable values with rcalizable frequency rcsolution mcasurements. 
The highest signal frequency of interest determincs the required sam- 
pling rate of the data. The sampling rate and the time window length 


determine the required dimension of the data input matrix. With these 


J; 





eonsiderations in mind, the following parameters were ehosen in digitiz- 
ing the data: 
‘Sampling Rate of 512 samples/seeond yields Nyquist frequeney 
ot Cooma: 
Input data low-pass filtered to 250Hz to prevent aliasing errors; 
Time reeord length of 20 seconds for FFT resolution of .05 
Hz / Bin; 
Hanning window funetion applied in time domain to reduee 
sidelobes associated with rectangular data windows. 
2. Power Speetrum 
With the above chosen data parameters, the input matrix to the 
POWER Subroutine* [3] was dimensioned to 20, 480 bins. This eor- 
responds to two bins per data point (1 cosine term and 1 sine term for 
10, 240 data points) or 10, 240 Complex Fourier eoeffieients. The real 
data values were loaded into the odd (cosine or real) bins and zero's 
were provided for the even (sine or imaginary) bins. On return from 
the FFT Subroutine [6] to the POWER Subroutine, the eosine term is 
squared and added to its assoeiated sine term squared thus representing 
the power at some diserete frequency. ‘This discrete frequency is an 
integer multiple of the FFT proeessing resolution (the reeiproca) of the 
time reeord length). The first bin, Y(1), contains the zero frequeney 
power eoeffieient, the second bin eontains the .05 Hz power eoefficient, 
and so on in .05 Hz steps. Thus bins У(1) through Y(5001) now contained 
the power eoefficients from zero to 250 Hz. 


“Subroutine STATS and POWER were written by Arfman, J. F., Jr., 
and were adapted for use in this program. 








3. Search Window and Background Noise Suppression 

With over five thousand power coefficients available for track- 
ing, it is convenient to limit the range of search to those that offer the 
best data for doppler tracking. A signal having favorable characteristics 
is first identified, i.e., frequency above 100 to 150 Hz, apparently 
stable on a one Hz resolution processor, and of sufficient strength to 
maintain a favorable E ratio for a mre of the time. 
Then a search window is designated by the adjustable parameters FL 
and SRBW, where FL is the lower frequency limit of the search window 
and SRBW is the search window bandwidth. Both parameters must be 
some integer multiple of the FFT resolution and would usually be kept 
at whole numbers for convenience. The total number of power coefficients 
in the search window would equal N=20(SRBW) and are designated YA(1) to 
SAN). 

Ideally when the frequency to be tracked (FTRK) has been iden- 
tified, FL and SRBW should be chosen such that FTRK is equal to FL 
plus (1/2)SRBW and SRBW is as large as is reasonable to keep the ratio 


of signal power in the window to the total power in the window much less 


L M 
than one. If Ж ҮА(К) << 2 YA(J) for L signal bins and M noise 
IT Aes 


bins with L4MzN, then the summation as I goes from one to N of YA(I)/N 
is a reasonable approximation of the ambient noisc level around FTRk [4], 
[5]. This calculation is performed and set equal to YNOIS. Then YNOIS, 
or some function of 1t dependent on the signal-to-noise ratio, is sub- 


tracted from cach power coefficient in the search window. Under the 
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assumption that the noise power spectrum is constant (white noisc) 
within SRBW, only the signal power spectrum remains in the window. 

4. Frequency Tracking to Follow the Doppler Shift 

The frequency selected to be tracked will be the strongest sig- 

nal in the search window when the first transform is computed. The 
bin number of that signal is saved as TBIN(JJ). In each succeeding 
transform the strongest signal (largest coefficient) in the search window 
is designated TBIN(JJ) and is compared to TBIN(JJ.1). A bin shift 
greater than 10 is determined to be a new or different discrete signal. 
lhe new power coefficients around the last center bin are then examined 
for discrete signal energy. The maximum value found would be con- 
sidered the signal energy being tracked if its bin number is equal to the 
last center bin number plus or minus 1, and its magnitude is greatcr 
than a chosen threshold minimum. In this problem this value was set 
to three times the ambient noise value. If the maximum coefficient 
found is more than one bin from the last center bin, the old center bin is 
examined. If its coefficient is larger than at least one of its sidelobe 
bin values and at least half as large as the maximum found and larger 
than the chosen threshold minimum value, the old signal bin is con- 
sidered to still contain the signal being tracked. If any of these tesi 
fail, the new maximum found is compared to the chosen minimum 
value. H the new maximum is larger, it is accepted as the signal being 
tracked; if not the contact is considered lost. Figure l possibly 


clarifies the above process. 
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Once the bin containing the largest portion of the signal energy 
being tracked has been established, the frequency of that signal is 
refined by computing the first moment of the three bin power "hump" 
around it. If the detected Signalis an integer ETE of the FFP 
resolution, adjacent bin spill over will be minimum and the refining 
calculation would make no correction to the frequency of the bit 
associated with the maximum power. When the detected signal is not 
an integer multiple of the FFT resolution the refining calculation adjusts 
the frequency to a mid-resolution value corresponding to the first mo- 
ment of the hump. See Figure 2, The frequency corresponding to bin 
number K is stored in the YL matrix as YL(K). The frequency refining 
equation used was 


FTRK(JJ)-YL(K)4RESOL*(YA(K-1)- YA(K- 1))/(YA(K41) 
SONS N (3) 


where 
FTRK=the discrete frequency being tracked, 
RESOL=the resolution of the FFT processing, 
YA(M)=the power coefficient of the Mth bin. 
Letting YL(K)=220. 25Hz and using the sample values from figure 2, 


| 


FTRK(JJ) 20ےے‎ Votan] ( 0 AE E 


i! 


ЕРЕ А لوا‎ (25) 
=: 220. 25017, - . 01252 
= 032 3 ٣ 


220, 2321, rounded) 
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Without this calculation the returned value of the frequency being tracked 
would have been 220.250Hz. A five bin computation was found to offer 
no significant improvement over the three bin computation due to the 


increased probability of noise or adjacent discrete signal influence. 


С. DESIGNATING PROBLEM VALUES 
In computing the sound propagation speed for the following examples, 
nominal values were used because the true values were not available. 
Values used were: 
0 


Temperature = 48 F 


Depth 60 ft 


Salinity рр 

(The sound emitter on the data tape appeared to make a slow turn 
into the buoy for the first 4 minutes of the tape and then maintain a 
constant course at a speed of approximately 4.5 knots. The closest- 
point-of-approach (CPA) to the hydrophone was approximately 300 to 700 
yards and 20 minutes of data was utilized with CPA at about the 11 to 12 
minute point on the tape. A 20 Hz search window bandwidth was used 
with several different lower frequency limits on the window. It was 
hound that one transform per minute comfortably contained the tracked 


frequencies within the 2] bin tracking window in all cases for a slow- 
q 


moving, near-field emitter. 
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DDATA OUTPUT 
Ia l (See Table |) 
With FL set to 40 Hz and SRBW at 20 Hz, the rest frequency 
of the tracked signal was approximately 50. 13 Hz. The total doppler 
shift was only 0.15 Hz and was too low to provide good tracking data. 
2. Run2 (See Table 2) 
For this run FL was set to 220 Hz with SRBW still at 20 Hz. 
Two strong discrete frequency signals were within the search window 
bandwidth, one's rest frequency at 233.87 Hz and the other at 223, 22Hz. 
The maximum power coefficient changed between these two signals 
several times. Inthe initial transform the tracker locked on to the 
lower of the two frequencies and a continuous track was maintained even 
though the signal faded out and then carne back. This run gave the fre. 
Buency tracker a good check-out in that (a) the maximum coefficient in 
the search window had considerable variations between discrete signals 
(b) a noise spike occurred in the tracker window at time=430 (FILE=220) 
and (c) the tracked signal faded completely out at timez1090 (FILE-550) 
and was regained one minute later. The total doppler shift was approx- 
imately 0.6 Hz and provided good data for doppler tracking. The computer 
output for this run is included in this report in its entirety. 
3. Run 3 (See Table 3) 
For this run FL was set to225 Hz and SRBW remained at 20 Hz. 


The only significant signal in the search window was the 233.87 Hz signal 
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Time Center Ratio: Frequency 


(seconds) Bin Signal/Noise Tracked 
10.0 204 155/66 50. 155 
20.0 205 1759/38 50.190 

130.0 205 1619/59 50.185 
190.0 205 2001/54 50.203 
250.0 205 1620/39 50.180 
310.0 205 1527/49 50.192 
370.0 204 878/40 50. 197 
430.0 205 247/42 1 
490.0 205 9320/78 50.178 
550.0 205 13280/68 50.191 
610.0 204 56746/254 50.168 
670.0 203 38878/144 50. 105 
730.0 203 4 ٦ 507055 
790.0 2035 47640 /227 50.081 
850.0 202 6941/55 50.060 
910.0 202 9654/71 50.063 
070.0 202 18841 /125 5605055 
1030.0 202 17558/111 507057 
1090.0 202 17546/104 50.054 
1150.0 202 14377/87 50.050 


Rest Frequency approximately 50. 13 Hz 


TABLET. 
11 1 ٢ TRACKER OUTPUT WITH FL= 50 Hz, SRBW=20 На. 
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Time Center Ratio: Frequency 


(seconds) Bin Signal/Noise Tracked | 
10.0 70 69/6 22 2 27 
70.0 71 112/5 22ے‎ | 

130.0 71 61/5 223.495 
190.0 71/2 52/5 223.556 
250.0 71 254/6 223.489 
310.0 70 218/6 225 ٦ 
370.0 69 490/9 223.420 
430.0 69 64/8 223.401 
490.0 70 743/24 223.434 
550.0 1! 2131/41 223.488 
610.0 69 10429/68 223.396 
670.0 62 4795/50 223.003 
730.0 60 2554/26 222.960 
790.0 61 2119/16 222. 995 
850.0 60 322/14 222.938 
910.0 59 245/10 LlLl, OBL 
970.0 59 60/4 222.898 

1030.0 59 112/6 222.886 

1090.0 60 12/5 Lost Track 

1150.0 59 115/5 222.901 


Rest Frequency approximately 223.22 Hz 


TABLE II. 
PRECQUENEY TP CIEELR OUTPUT WITH FL=220 Hz, SRBW=20 Hz. 
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Time Genter Ratio: Frequency 


(seconds) Bin Signal/Noise Tracked 
10:0 183 61/4 234.100 
70.0 184 42/4 234,151 

150.0 185 21/4 2342193 
190.0 185 77/5 234.224 
250.0 184 167/5 234.169 
310.0 184 160/5 23534 ۶ 
о 0 183 829/1 234. 102 
430.0 184 708/8 234. 144 
490. 0 184 3295/19 2947132 
550.0 184 4824/32 2341910) 
610.0 183 4669/30 234.079 
670.0 176 6049/29 233.756 
750.0 174 906/13 233 4 
790.0 174 766/7 283 A 
850.0 172 1082/11 233.563 
910.0 172 1020/8 2550559 
970.0 172 LIA 2290 9 

1030. 0 172 357/5 233,545 

1090.0 172 135 /5 233 30 

1150.0 1 23.05 233 525 


Rest Frequency approximately 233.87 Hz 


A JII. 
J P OU C Y TI CEERP OUEEUT WITII FIS 225 Hz, SRBW-20 Hz. 
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and this was tracked with only à minimal amount of decision require. 
ments around CPA. The total doppler shift was approximately 0.65 Hz 


and provided good data for doppler tracking. 


ШІ TARGET TRACKING 

Since no bearing information was available corresponding to the 
frequency data, the IBM 360 Computer was used to generate simulated 
real time bearings. This was accomplished by programming a pseudo 
"true track'' into the computer so that the target speed and the CPA time 
were matched as wellas possible. The initial starting position was 
X=500 yards, Y=-1050 yards at Fie T(1)=190 seconds. The heading 
and speed used were 90 degrees and 4.5 knots where zero degrees 
corresponds to the positive X axis. The computer then computed "true 
bearings'' to the target and added random noise values of a programmed 
standard deviation, SA, to these values. Thus the true track plotted 
by the computer is only a best estimate of the targets relative position 
from the buoy. 

The X- Y Filter developed by Mitschang [1] was used to process the 
data and develop the target track. To initialize the filter and commence 
tracking, the measured bearing must have changed by at least three 
standard deviations of the measurement noise and the doppler shifted 
frequency must have decreascd from lts valut at the start of the initializa - 
tion process. A modification to the filter was attempted so that while 


initializing, an increasing frequency would be detected as a mancuvcering 





target. Time, bearing, and frequency values occurring with the maxi- 
mum ‘frequency detected would then be retained as the initial values to 
be used when the bearing shift and non-increasing frequency criteria 
were again satisfied. This would delay the start of the tracking process 
but would prevent data from a maneuvering target being used to initialize 
the state equations based on a constant heading, constant speed target. 
However, indexing problems were encountered and ven not resolved 
in time to include the modification in this report. 

The following difference equations developed in Chapter 3 of [1] 


were used to compute the initialrange, heading, and rest frequency of 


the target. 
а (Мао о DELE) 0)8 ۰“ i ) (4) 
0si-0ave-180 - arcsin((-VP- DELF)/(f(last): DELO- Vsi)) (5) 
Е0і-Ғаче (1 - (уві/УРусов(Өві - баче)| (6) 
where 
ҮР =speed of sound propagation 


DELF =doppler shift of the frequency 


DELT =time required to meet initialization criteria 
f(last) =last frequency measured for en 
D ELO -change in bearing from sensor to targct 

Osi =inilial estimate of target heading 

Dave =avcrage bearing for the initialization period 
Vsi =initial estimate of target speed 
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ЕО: =initial estimate of the rest frequency 
Fave =average frequency for the initialization period 
Then based ona constant heading, constant speed target, the initial 


filter state equations are: 


X(1) =X position = R cos 6 (7) 
X(2) =X velocity = Vsi cos 651 (8) 
X(3) =Y position = R sin 03۶٤٥ (9) 
X (4) = У velocity = Vsi sin Qsi (10) 
Х(5) =Rest frequency = FOi (11) 


The initial covariance matrix values are computed by the "direct 
method"! developed in Chapter 4 of [1] and an extended Kalman filtering 
technique = used to refine the estimates of the target's position, 
heading, speed, and rest frequency as the problem progresses. 

Since the target was maneuvering for the first four minutes of the 
tape, initialization of the target tracker was commenced at time 
T(4)=190 seconds, with the values from Table 3 as input data. The 
average initialization time, AVKJ, varied with the standard deviation of 
the bearing measurement noise, SA. This time varied from about 4 


minutes for SA=5 degrees to about 8 minutes for 5٨۵۸-30 
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hie oligo OBTAINED FROM REAL FREQUENCY INPUTS 


Onee the frequeney tracker has loeked on to a stable but doppler 
shifted frequeney, the inputs required of the operator are an estimate of 
the rest frequeney, FO, and an estimate of the target's speed, AVSK. 
These.are ''priming' values used in eomputing the initial heading of the 
target and affect the ultimate aeeuracy of the eomputed traek. A third 
variable affeeting the initialization is the standard deviation of the bear- 
ing error, SA. As SA increases, the filter takes longer to start traeking 
and has more data on whieh to base its original estimates. 

The program was run for various eombinations of the three variables 
FO, AVSK, and SA. The true value for FO was approximately 233.67 112 
and for AVSK was 4.5 knots. An estimated FO greater than the deteeted 
frequeney causes a 'down-doppler'' solution for the initial eourse. Con- 
versely an estimate of FO less than the detected frequeney results in an 
"up-doppler" solution for the initial course. Eaeh combination of var- 
iables was looped so that one hundred runs were made with that set of 
inputs and the averaged filtered positions were plotted. The "+" symbol 
on eaeh plot represents the eomputed position of the target at time T(1) 
plus one-half of the initialization period. The range is eomputed by 
Equation (4) and the bearing is the average bearing for the initialization 


period. The true points from whieh the noisy bearing measurements 
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were made are plotted with the symbol ''*"', Averaged filter positions 
are plotted with an "X€". The ascending order of precedence on each 
plot is +, *, X. Figures 3 through 24 show the performance of the system 
with real, high resolution frequencies as inputs. F(DET) signifies 


initial frequency detected. 
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"Y, ZE@NCLUSIONS 


The automatic frequency tracker developed in this report has been 
shown to be adequate for tracking a doppler-shifting frequency ina 
favorable signal-to-noise ratio environment. The width of the tracking 
^window-.-would:havwe.to.he.increased for a higher speed target unless the 
transforms were computed at closer time intervals. A loss of signal 
would result in the last center bin frequency being retained as the signal 
being tracked unless a spurious signal with more power than the chosen 
minimum appeared in the window. Therefore, under marginal signal- 
to-noise ratio conditions, an operator would be required to monitor the 


'' condition 


frequency tracker output to insure that after a "lost track 
occurred, any "regain track'' operation was indeed occurring on the 
previous signal being tracked. 

The procedure of using Equation (3) to extend the frequency resolu- 
tion of the output above that of the FFT processing is dependent upon 
having a reasonable approximation of the ambient noise level. Comput- 
ing a new estimate in each transform results in higher power levels 
than the true ambient noise level being used when the target is near CPA. 
Since the ratio of two numbers, A/B, is not equal to the ratio (A-K)/ 
(B-K) where K is a constant, the subtraction of the corrupted estimate 


of ambient noise level causes an error in the value of the resolved 


frequency. However, in comparing the amplitudes of the signal and 


Бе 





noise estimates around CPA (from Tables I, II, and III), the signal is 
150 to 200 times as large as the corrupted ambient noise estimate. The 
resultant error on the resolved frequency is on the order of magnitude of 
1074^Hz and was considered negligible. As the range at CPA increases, 
the signal power level decreases and its corrupting influence on the 
ambient noise level approximation is decreased. Thus as long as the 
increased average power level in the window is caused bs the target 
being tracked, no appreciable error occurs and the resolution of the 
overall system has been extended beyond the resolution of the FFT 
processing. 

With real frequency inputs and reasonable initial estimates of the 
rest frequency and target speed, the Mitschang Filter provided target 
tracks in excellent agreement with the known emitter locations and 
probable emitter locations. Thus a true single sensor fix appears 
possible using existing hardware and techniques. Further testing is 
needed where exact target locations are known and real measured bear. 


ings are available corresponding to the real frequency measurements. 
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